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We present highlights of the spectrum of mesons and baryons calculated using NRQCD for heavy quarks and 
tadpole improved clover action for the light quarks. 
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1. HEAVY-LIGHT MESONS 

The results for heavy-light mesons are sum- 
marized in Table |lj. The results presented here 
are based on the same statistical sample, consist- 
ing of 102 quenched configurations at f3 = 6.0 
with lattice size 16 3 x 48, as used in our study of 
Jb and f Bs §. The details of the NRQCD ac- 
tion, the evolution equation for calculating the 
heavy quark propagator, the method used for 
setting the lattice scale, and the fixing of light 
and strange quark masses are given in JD. The 
lattice scale, determined from M p , is 1.92(7) 
GeV, and we consider the range 1.8 — 2.0 GeV 
to determine the associated uncertainty. The 
bare b quark mass used in the heavy-light anal- 
yses is aM® — 2.31(12). From this we esti- 
mate m^im^S) = 4.35(10) (-^) (10) GeV. For 
hadrons containing strange quarks, the central 
values correspond to fixing m s using Mk, and 
for errors we use the difference between it and 
using itl s (Mk*)- Details of the analyses of the 
spectrum will appear soon in |2| . 

Overall, our estimates are in rough agreement 
with experimental data, the one exception being 
the hyperfine splittings as discussed below. To 
understand the various mass splittings, we use the 
following qualitative picture in which the mass of 
a heavy-light hadron is considered to be a sum of: 

• the pole mass of the heavy quark (Mh) 
which is ~ 1.5 GeV for the c quark and 
- 5.0 GeV for the b; 



state (n J 1 ") 


Lattice 


Expt. 


heavy-light mesons 


B 1(0-) 

2(0-) 

B* l(l-) 

B* 1(0+) 
B 1 1(1+) 
Bl 1(2+) 


5296(04) (+3) 
5895(1 16) (+J!j) 
5319(02)(+°) 
5670(37)(±J2) 

5726(38)(+J°) 
5822(45) (til) 


5279 

5860(*) 

5325(1) 

5698(12) 

5779(*)i 


heavy-strange mesons 


B s 1(0-) 

2(0-) 
B* s 1(1") 

B* s0 1(0+) 
B sl 1(1+) 

B* S 2 1(2+) 


5385(15)(;?)( + -) 
5935(57)(tl)(i 9 ) 
5412(14)(^)(±g°) 

5742(27)(+i 4 )(+^ 5 ) 
5804(31)(+^)(±J 6 ) 
5878(26)(±1)(± 1 ) 


5375(6) 
5422(6) 
5853(15) 



*This work was done in collaboration with A. Ali Khan, 
T. Bhattacharya, S. Collins, C. Davies, C. Morningstar, 
J. Shigemitsu and J. Sloan, and supported by the DoE 
Grand Challenges award at the ACL at Los Alamos. 



Table 1 

Mass estimates in MeV for various meson states. 
The b quark mass is fixed using the spin-averaged 
B(1S). The first error of the lattice data is sta- 
tistical (this bootstrap estimate includes uncer- 
tainties due to extrapolations in quark masses), 
the second represents the scale uncertainty due 
to varying a^ 1 between 1.8 and 2.0 GeV, and for 
the strange mesons we quote a third error asso- 
ciated with the uncertainty in fixing the strange 
quark mass. Preliminary experimental values are 
denoted by asterisks. The experimental numbers 
quoted against the 1+ states correspond to broad 
resonances that are possibly a mixture of the dif- 
ferent P states. The corresponding lattice esti- 
mates are an unresolved combination of the 1+ 
and 1+' states. 



• the constituent mass of the light quarks 
which is approximately 300 MeV for the u, d 
and 450 MeV for the s quark as inferred 
from the octet and decuplet light baryons. 

• for orbitally and radially excited states, an 
excitation energy of the light quark, which 
we expect to be of the order of Aqcd ', 

• the 0(AQ CD /Mh) contributions are the ki- 
netic energy of the heavy quark and the 
heavy- light hyperfine energy E a - H . a - l rs 46 
MeV, inferred from the experimental B* — B 
splitting; 

• and a residual binding energy E be encapsu- 
lating the remaining interactions which we 
expect to be small (0(Aq CD /M^)). 

To isolate individual terms and estimate their size 
and dependence on the quark masses we construct 
different linear combinations of meson and baryon 
masses. 

The spin-averaged splitting M b b — M B d should 
be dominated by the difference of the strange and 
light quark masses. We find 90(9) [±\){tf) MeV 
(the first error includes statistical and extrapo- 
lation in quark masses, the second is due to the 
scale uncertainty, and the third is the variation 
if m s is set using Mr-* instead of Mr-) whereas 
the experimental value is 96(6) MeV. The split- 
ting shows no significant dependence on the heavy 
quark mass, consistent with the expectation that 
the change in the difference in the kinetic energy 
of the heavy quark is small. 

By a judicious combination of operators and 
sources for quark propagators we are able to ob- 
tain a signal for P-wave states. The estimates 
for 3 P2 and 3 Pq are unambiguous as the op- 
erators used to probe these states do not mix 
with other states. Thus, their masses and the 
splittings B% - B% = 155(32) (± 9 13 ) MeV and 

B* s2 - B* s0 = 136(23)(±i§)(±2) MeV, are predic- 
tions since the P states have not been resolved 
experimentally. We have used 3 Pi and l P\ op- 
erators in the LS coupling scheme to probe the 
physical 1 + and 1 +/ states. In this case each oper- 
ator has a non-zero overlap with the two physical 
states; consequently to get the physical masses it 



is essential to get a signal in both the direct and 
mixed correlators. Since the data show no reliable 
signal in the mixed correlators, we do not have 
predictions for the masses of the physical states. 
In Table [y, we quote the result from the 3 P\ cor- 
relator as a rough estimate for the 1 + state. 

Radial and orbital splittings are expected to 
be dominated by the difference in kinetic en- 
ergies of the heavy and light quarks. Of 
these, the light quark contributes the most 

°(^QCD I ™ constituent) ~ 0(Aq C d)- We find 

602(86)(j|jj) and 559 (55) (±j&) (±$ 2 ) MeV for the 

2 1 So — 1 * So splitting in the B and B s systems 
respectively. The preliminary experimental value 
for the B is 581 MeV pj. For the spin-averaged 
IP - 15 splitting we find 457(31) (±35) MeV for 
the B, and 428(27) (i^X^) MeV for the B s . 

The results for the hyperfine splittings are 
not as encouraging. We find AE(B^ — 
B d ) = 24(5)(±|) MeV and AE(B* S - B s ) = 
27(3)(±|)(±o) MeV, i.e. roughly half the exper- 
imental values, 46 and 47 MeV respectively. Fur- 
ther work is required to clarify whether this dis- 
crepancy is due to the quenched approximation 
or due to an underestimate of the a ■ B term (the 
clover term) in the action. 



2. BARYONS 

Our results for baryons are summarized in Ta- 
ble y. The first splitting we comment on is 
M\ h — (Mjj + 3Mjj») /4. In this there is no con- 
tribution from the heavy quark mass and hyper- 
fine interaction E a - H . a - l , so it should be dominated 
by mass of the extra light quark in A. The exper- 
imental values, 311(10) and 310(2) MeV for the 
b and c systems respectively, support this. Our 
lattice estimates are A b -B = 370(67)(±2o) MeV 
and E b -B~ s = 392(50)(+o 5 ) MeV. Also, the data 
show little dependence on the heavy quark mass. 

In our picture the hyperfine interaction be- 
tween the light quarks should dominate the split- 
ting (2E/ l +4E£)/6-A/ l . A calculation of the a v cj[ 
term, assuming a simple non-relativistic model, 
suggests that this splitting should be 2/3 of the 
Delta-Nucleon splitting which is 293 MeV. The 
experimental results, (2S C + 4S*)/6 — A e = 212 



baryon 


expt. 


Our results 


A-likc {si = 0, j= 1/2) 


A h {udb) 
E h {Isb) 


5.624(9) 


5.679(71)CT 

5.795(53)(±? 5 )(± 5 ) 


E-like {si = 1, j= 1/2) 


E ft {Ub) 
E' h {Isb) 
fl h {ssb) 


5.797(8) 


5.887(49) ( + ^) 
5.968(39)(+ 2 3 °)(^ 4 ) 
6.048(33)(^)(± 34 ) 


E-like (si = 1, j= 3/2) 


E* {lib) 
E* h {Isb) 


5.853(8) || 


5.909(47) (+^) 
5.989(39)(ig)(^ 5 ) 
6.069(34)(t 3 «)(^ 5 ) 





a' 1 MeV 


M°a 


M 6 MS GcV 


2 3 S ± - l 3 5i 
2 So — 1 Sq 

^i - 3 s 1 


2313(99) 
2413(83) 
2424(133) 


1.76(7) 
1.69(5) 
1.69(8) 


4.07(24) 
4.08(19) 
4.10(30) 



Table 2 

Summary of b baryons masses in GeV (I stands 

for auorrf quark). 



MeV and (2E 6 + 4E*)/6- A b = 210 MeV (prelim- 
inary) are consistent with this and suggest negli- 
gible dependence on the heavy quark mass. We 
find 221(71) (i^) MeV at M b and no significant 
dependence on the heavy quark mass. The value 
decreases to 186(51) (+lf)(i? ) MeV on replacing 
d with s, in qualitative agreement with the exper- 
imental results in the charmed sector. 

The E ; * — Eft splitting is expected to be pro- 
portional to (1/M/j) as it should depend only on 
the heavy-light hyperfine interaction E ah . ai . Us- 
ing this heavy quark scaling one expects a value 
(B* - B)(E* - E C )/(P>* - D) « 46 * 66/140 « 22 
MeV. We find 19(7) (if) MeV, whereas the pre- 
liminary experimental value is 56(8) MeV ||. 
However, the experimental identification of the 
states is still under debate. 



3. HEAVY-HEAVY MESONS ONIA 

The bottomonia spectrum has been success- 
fully used for an independent determination of 
the lattice scale and m\, || . Here we present pre- 
liminary results based on the same data set as 
for the heavy-light analyses presented above. To 
extract the lattice scale we assume that the split- 
tings (Mr- -M r ), (Afy -M„J, and (M % -Mr), 
are independent of the heavy quark mass. We in- 
clude (M„> — M Vb ) as it has the best signal and 



Table 3 

Estimates of lattice scale and b quark mass from 

T. The final column gives Mf !s (Mf is ). 



we assume it is equal to (Mr' — Mr) on basis 
of the approximate equality in the charm system. 
Having fixed 1/a we determine M® by fixing Mr 
using the dispersion relation. These results are 
summarized in table [3| 

Typically, the results for mass-splittings with 
best control over statistical errors are obtained 
from fits to ratios of correlators. In our data we 
find marginally better signal from fits to individ- 
ual as compared to ratio of correlators (a con- 
sequence of chosing a smearing in heavy quark 
propagator generation that optimizes signal for 
heavy-light states). For the hyperfine splitting 
Mr - M m we find 34.2(3.6) MeV (ratio fits) 
and 30.4(3.2) MeV (difference of mass fits). If 
one scales the J/tp — r\ c splitting, 127 MeV, by 
Mj/^/Mr then one expects ~ 42 MeV. Thus, 
once again the hyperfine splitting is underesti- 
mated though not by as much as for Mb- — Mb- 

The determination of P state splittings has 
larger errors. We find l P\ — 3 P = 0(8) where 
3 P = (fP + 3 * 3 Pi + 3 * 3 P2(T) + 2 * 3 P 2 (P))/9; 
3 Pp- 3 P = -42(23); 3 Pi- 3 P = 13(23); 3 P 2 {T)- 
3 P = 28(14); and 3 P (E) - 3 P = -40(23) 
MeV. It is disturbing to find the T and E cu- 
bic representations of 3 P2 give such different es- 
timates. The corresponding experimental values 
are -40.2, -8.3, and 13 MeV. 
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